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Abstract 

The static vortex solution in Abelian Higgs model with small ratio of 
vector and Higgs particle masses is considered. Several formulae approxi- 
mating this solution are discussed. The accuracy of these approximations 
is tested by numerical computations. 
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1 Introduction 



Nowadays vortex solutions are found to be interesting in many areas of 
physics. Their investigations can help in better understanding of some phenom- 
ena in field theory, cosmology and condensed matter physics jjj . It is rather dif- 
ficult to obtain vortex solutions since one has usually to solve highly non-linear, 
very complicated differential equations. Therefore appropriate analytical and 
numerical methods must be worked out and applied to make some progress in 
this area. 

One of the simplest systems possesing vortex solution is the Abelian Higgs 
model. An exact vortex solution, in the form of infinite convergent series, was 
found for this model in the so-called Bogomolny limit (the masses of the scalar 
and vector bosons are the same) Q . In this limit the equations of motion reduce 
to the first order differential equations and the underlying methods can not be 
simply applied in general case. 

Recently in the Abelian Higgs model a polynomial approximation was ex- 
tensively used to investigate excited vortex Q , Q • This method gives analytical 
formulae simple enough to be applied in further computations. However it is 
necessary to estimate the error of this approximation and this can be done by 
a comparison with the numerical solution. 

The goal of the present paper is twofold. First we would like to test some ap- 
proximating formulae for the static vortex solution proposed in paper || . These 
analytic formulae were obtained in the limit n — > where kappa is the ratio of 
vector and Higgs particle masses. Our investigations are also limited to small 
values of k. Second we would like to present more precise analytical approx- 
imations. The accurate analytical and numerical results for the static vortex 
solution can be useful in many problems apart from the excited vortex. For 
example, equations governing evolution of a curved vortex (5| involve constants 
which are determined by the static vortex solution. 

Our paper is organized as follows. In Sec. 2 we shortly review the simple 
analytic formulae approximating the static vortex solution. As mentioned above 
these formulae were introduced and discussed in more detail in ||. We also 
compare the results obtained this way with numerical approximations to test 
the accuracy of the algorithm. In Sec. 3 we propose several improvements of 
this method to get more precise analytical approximations and numerical results 
valid in wider region than the previous ones. Finally in Sec. 4 we present some 
general remarks and conclusions summarizing our paper. 
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2 Simple approximations 



The Abelian Higgs model in 3 + 1-dimcnsional Minkowski space-time is de- 
scribed by the following Euler-Lagrange equations 

(d M + iqA^ + iqA^ + - ^) = 0, (1) 

d^ v = iq($*d u $ - - 2q 2 A»\<P\ 2 . (2) 

The Higgs field <&(x) is complex valued, the star denotes the complex con- 
jugation. The signature of metric of the space-time is ( — h ++). The static 
Abrikosov-Nielsen-Olesen solution [|| represents a straightlinear, infinite vortex. 
The unit topological charge vortex lying along the z-axis can be obtained from 
the equations by imposing on them the axially symmetric Ansatz 



A = 0, A 3 = 0, 

A 1 =sin(g) x(p) ~ 1 , (3) 

qp 

A 2 = -cos(6) x{p) -\ 

qp 



Here p = y/ (x 1 ) 2 + (x 2 ) 2 and 9 = arctan(x 2 /a; 1 ) are the standard polar coordi- 
nates in the (x 1 ^ 2 ) plane. 

The Ansatz (||) together with the rescaling of the p variable 

r = 2m 2 p, (4) 
leads to the second order differential equations for the fields F(r) and x( r ) 

F" + — -X-F+Uf-F 3 )=0, (5) 
r r l 2 

x"-^-k 2 F 2 X = 0, (6) 
r 



where prime denotes differentiation with respect to r and k = y is the only 
remaining free parameter. 

To guarantee the vortex solution to be regular on x 3 -axis and to have finite 
energy per unit of length in x 3 direction the above equations must be supple- 
mented by the boundary conditions 

F(0) = 0,F(oo) = l, (7) 
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X(0) = l,x(oo) = 0. (8) 

The asymptotic form of \ can be easily obtained from eqs. @,(^) simplified 
by putting F = 1. Thus we get 

/ 7TT 

Xasym{r) = c rKi(nr) w c J — exp(-Kx), (9) 



where K\ is the modified Bessel function 



and co is a constant. The asymptotic 



behaviour of F(r) strongly depends on k [3|. For n < ^ it is determined mainly 
by the term y^Fjr 2 and can be obtained from eqs. (||),(@) 

2 

Fasvml = 1 - k(1 C f An2y exp(-2 K r), (10) 

while for k > i the Higgs term (F — F 3 )/2 is more important. In this case eq. 
(^|) linearized in (F — 1) gives the following result 

/ 7T 

Fosj/ma = l + CiiifoCr) « 1 + ciW — exp(-r). (11) 

Here Kq denotes the zero-order modified Bessel function (7) and c\ is a constant. 

The asymptotic forms of F and x were used in || to get the approximate 
static vortex solution. In the neighbourhood of r — the fields F and x were 
approximated by low order polynomials obtained as power series solutions of 
eqs.(||),(^|). These polynomials were smoothly matched with the appropriate 
asymptotics at some point r — ro i.e. the functions F(r) and x( r ) were required 
to be continuous at r — tq together with their first derivatives 

F po iy(r ) = F asym (r ), F' poly {r ) = F' asym {r ) (12) 



Xpoly(ro) = Xasym(ro), x' po ly( r o) = xLymM- (13) 

In the simplest version proposed in the static vortex solution was approx- 
imated by the following formulae 

F(r) = { {^-y^ * r<r ° (14) 
w [ 1 if r > ro, 

r s _ / 1 - hx2r 2 + i X4 r 4 - ^xer 6 if r < r 
X[T> ~ { corK^r) if r > r , [L *> 

involving four constants r , fx, X2, c o- These constants were fixed by applying 
the matching conditions described above. The other constants are given by 
recurrence relations (jl^) below. 
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In order to get more accurate solutions the formula for F in the region r > tq 
was replaced with a more subtle one 



F= r w' (16) 

obtained from eq. (||) simplified by neglecting the terms with the derivatives of 
F. In this case the polynomial approximation of F must be completed with the 
term proportional to r 5 

F = fir-^hr 3 + ±f 5 r 5 : (17) 

while the formulae for x remain unchanged although the values of the particular 
parameters are different. The equations (d),© l ea d to the following recurrence 
relations for coefficients of the polynomials 



h = \{\ + X2)h + \{\xi + \xl + fl)h, (18) 

X6 = 5 K 2 / 1 (2/ 3 + 3 X2 /i). 

In Fig. (la)-(lc) we have compared the described above approximation 
of the Higgs field F with its numerical values obtained by applying standard 
algorithms for stiff differential equations ||. As we are interested mainly in 
small values of k we have limited ourselves to K = 0.02, k = 0.1 and n = 0.2. 
Since the values of the field x obtained from the numerical computations and 
approximate formula differ very slightly we have plotted their differences ( Fig. 
(2) ) and the numerical values themselves ( Fig. (3) ). The numerical values of 
the free parameters /1, X2, tq, cq are given below. 



(F aS y m 



K 


h 


h 2 


ro 


co 


0.02 


0.6505427 


0.00157722 


2.305767 


0.020004 


0.1 


0.6646855 


0.02362296 


2.256706 


0.1005027 


0.2 


0.6929167 


0.06904716 


2.164762 


0.2040623 



(Fasyrn = V 1 ~ 2(xA) 2 ) 



K 


fi 


h 2 


ro 


CO 


0.02 


0.4285536 


0.001431718 


2.106883 


0.0200066 


0.1 


0.443162 


0.02028046 


2.068759 


0.1007802 


0.2 


0.469365 


0.05754442 


2.000099 


0.2056864 
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Let us note that the approximate formula for the function x is quite good 
while the approximation of F is much worse. This is the price for simplicity 
of the analytical expressions. The field F tends very quickly to its asymptotic 
form and such behaviour can be hardly described by simple analytical formula. 



3 Improved approximate solutions 



The main defect of the approximate formulae considered in the previous sec- 
tion is the behaviour of the held F in the region of intermediate values of r 
particularly in the neighbourhood of the matching point ro. One can try to 
improve that approximation by using the higher order polynomial solution for 
the functions F and x m t ne interval (0, ro)- However the practical effect of 
such improvement seems to be rather small. A better accuracy can be reached 
by changing the approximation for the function F in the region (ro, oo). The 
more accurate asymptotics [|J is given by 



F asym = sll-2{^} + Cl Xo(r). (19) 

This formula involves a new parameter c\ and an extra condition is necces- 
sary to determine it. Therefore we have used an additional matching condition 
ensuring the continuity of the second order derivative of F at r = ?*o 

F poly( r o) = K'sym( r 0)- (20) 

We were able to satisfy the matching conditions if the polynomials approx- 
imating the functions F and % were of order fifteen and fourteen or nineteen 
and eighteen, respectively. We do not present them as their forms are very 
complicated and the numerical results not excellent as is shown in Fig. (4) and 
(5). 

There is also another simple posssibility to determine the values of five pa- 
rameters Co, C\, /i, h%, ro. One can solve the four matching conditions (|l2|) , ([l3|) 
for the fixed value of the radius ro and repeat this procedure for several values 
of ro in some interval. Thus we have the four parameters Co,Ci,/i,/i2 as the 
numerical functions of ro . The last step was to compare the approximations 
for F an x obtained this way with numerical calculations and fix the value of 
ro which gives the best fitting. It turned out that in this case it was enough 
to approximate F and x by polynomials of order five and four respectively, see 
eqs. (IIeI), ([l?]), (18) with neglected xe ~ terms. The differences between ap- 



proximated and numerical values of F and x are presented in Fig. (6) and (7). 
The numerical values of the parameters /i,X2, Co, ci are given in the following 
table 



G 



K 


h 


hi 


CO 


Cl 


0.02 


0.3811035 


0.001427667 


0.0200068 


1.368773 


0.1 


0.3885089 


0.02025245 


0.1007819 


1.530144 


0.2 


0.3998954 


0.05794554 


0.2052965 


1.851151 



while ro = 2.5. The vortex solutions considered so far were obtained by smooth 
matching of some polynomials approximating the vortex core with asymptotic 
formulae valid in the outer region. More accurate approximations can be ob- 
tained by dividing the whole area into more pieces and approximating the solu- 
tion in each sector separately. We have chosen k+1 points < r < r\ < ... < r k . 
The central part of the vortex in the interval (0, r ) was approximated by 

1 1 (—l) n 

F = frr - 3^ + -*r« + ... + (21) 

1 1 (—'\Y n 

X = 1 - ^ 2 + ij X4. 4 + + \^TX 2n r 2 -, (22) 

while in (rj , rj+i) for j = 0, 1, 2, ..k — 1 we have used the truncated Taylor series 
expansions 

F = f 0j + fij(r - rj ) + i/ 2j (r - r 3 f + ... + I/ ni (r - r^, (23) 



X = Xoj + Xij(r ~ r s ) + j { X2j(r - r 3 ) 2 + ... + ^Xnj{r - r^. (24) 

In the region (rfe,oo) the previous asymptotic formulae (|^), (19) were applied. 
We have required the functions F and x together with their hrst derivatives to 
be continuous in the matching points ro, r\, r k . These conditions and eqs. 

arc enough to determine all the coefficients fj, Xh fih Xij > °0t °i- Let 
us note that passing from r = to r = r k resembles the process of analytic 
continuation and is rather easy to perform. The main difficulty is to bind these 
solutions with their appropriate asymptotics. 

The above formulae were applied in two ways. In equations (^3|), ( |2~i| ) we 
have firstly put the expansion order to 4 (n = 2, rii = 4) and have divided the 
whole region of r into three pieces (r = 2,r\ = 3). In the region (0,r ) the 
formulae (pl|),(p2]) for n = 2 reduce to (|T5|), (|l7|), ( |l^ ) with neglected x 6 — terms 
while in the interval (ro,ri) the following recurrence relations are valid in eqs. 

I.© 

7 —fa . fo h . foXo 

~ 272- , XX 
X2 = « f X0 + , 
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fi 2/qxoXi XoA _h_[± + 3 /o/i _ 2 /oXo 
r§ r§ r 2 2 rg ' 



X3 = J + K faXi + 2k fa f iXo H , 

r o r o 

7 2/i 6/oXo /3 /2 . 3/ 2 /o , 0T27 

2/0X0X2 , 2/oXi , /2X0 , 4/iXoXi , 2/2 
8/0X0X1 4/iXo 



'O 'O 'O 'O 'O 



r o r l 



~ 2x2 . 2xi 272- , ,i 2?7~ 

X4 = 2- + — + « /oX2 + 4/t^/o/lXl 

r r 

+ 2 K 2 ./ xo/2 + 2k 2 xo/? + — , 

ro 

where the second index (j = 0) is omited for simplicity. The numerical values 
of the remaining parameters are found from the matching conditions and they 
are given in the following two tables 



K 


fi 


h 2 


co 


Cl 


0.02 


0.3980053 


0.1392637 


0.02000953 


1.273493 


0.1 


0.4078238 


0.01947355 


0.1011135 


1.34181 


0.2 


0.4238514 


0.05537332 


0.2079306 


1.479463 



K 


fo 


fi 


Xo 


Xi 


0.02 


0.6551700 


0.2454606 


0.9973415 


-0.002531821 


0.1 


0.6676333 


0.2496415 


0.9643793 


-0.03229428 


0.2 


0.6863076 


0.2557591 


0.9036254 


-0.08200263 



The rather simple but quite accurate analytical approximation of the static 
vortex solution was obtained this way, as is presented in Fig. (8) and (9). 

In the second case we have tried to get the possibly most accurate numerical 
results. Therefore we have put large no = 9, nj = 10 and very small value of 
r J+1 — rj = 0.01 for r = 1,77- = 20 . The corresponding results are shown in 
Fig. (10) and (11). 



4 Remarks 

In the present paper we have considered several formulae approximating the 
vortex solution in the Abelian Higgs model. We started with simple analytical 
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formulae presented in H and compared them with numerical computations. It 
turned out that the approximation of the Higgs field in the neighbourhood of 
the vortex core is rather rough and should be improved to get more accurate 
results. We have tried several methods to reach this goal. 

First of all we have changed the formula describing the Higgs field in the 
outer region. This formula should not be interpreted as the better asymptotics 
only. Perhaps more important is the fact that this expression involves a new 
free parameter which let us improve the whole algorithm. We have used this 
possibility in several ways. At this point it is worth noting that the relative 
error of the function F — 1 with F given by the approximate formula (jlq) does 
not tend to zero for large r. This can be easily seen by comparing eqs. (||), ( |To| ) 
and @. 

Our first trial to improve the accuracy of the approximation was the al- 
gorithm with an additional matching condition ensuring the continuity of the 
second derivative of the Higgs field. The next possibility we have tried was to 
solve the matching conditions in some fixed point and repeat this step several 
times in different points to choose finally the best matching point on the basis 
of numerical results. 

At last we have modified the algorithm by solving the equations of motion 
approximatelly as the truncations of the Taylor series expansions around an 
arbitrary point. We have used these solutions in the manner resembling the 
process of analytical continuation. This way we have obtained both: our best 
numerical approximations of the static vortex solutions and quite simple but ac- 
curate analytical formulae generalizing those from Q. It was possible because 
the final version of the algorithm turned out to be very flexible and could be ap- 
plied to reach apparently different purposes: analytical simplicity of expressions 
and numerical accuracy of computer calculations. 
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Figures 



la. The approximate and numerical values of the field F,k = 0.05. 
lb. The approximate and numerical values of the field F, k = 0.1. 
lc. The approximate and numerical values of the field F,n — 0.2. 

2. The differences between approximate and numerical values of the 
field X for k = 0.05,0.1,0.2. 

3. The numerical values of the field x for k = 0.05, 0.1, 0.2. 

4. The differences between approximate and numerical values of the 
field F for k = 0.05, 0.1, 0.2 (second order derivative method). 

5. The differences between approximate and numerical values of the 
field x for k = 0.05,0.1,0.2 (second order derivative method). 

6. The differences between approximate and numerical values of the 
field F for n — 0.05, 0.1, 0.2 (fixed matching point algorithm). 

7. The differences between approximate and numerical values of the 
field x for k = 0.05,0.1,0.2 (fixed matching point algorithm). 

8. The differences between approximate and numerical values of the 
field F for k = 0.05,0.1,0.2 (two matching points method). 

9. The differences between approximate and numerical values of the 
field x for k = 0.05,0.1,0.2 (two matching points method). 

10. The differences between approximate and numerical values of the 
field F for n — 0.05, 0.1, 0.2 (analytic continuation alghorithm). 

11. The differences between approximate and numerical values of the 
field x for k = 0.05, 0.1, 0.2 (analytic continuation alghorithm). 
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Fig. 1a 




Fig. 1b 




Fig. 1c 
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Fig. 5 
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Fig. 8 
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